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We obtain the band structure of a particle moving in a magnetic spin texture, classified by its 
chirality and structure factor, in the presence of spin-orbit coupling. This rich interplay leads to a 
variety of novel topological phases characterized by the Berry curvature and their associated Chern 
numbers. We suggest methods of experimentally exploring these topological phases by Hall drift 
measurements of the Chern number and Berry phase interferometry to map the Berry curvature. 


I. INTRODUCTION 

Strong spin-orbit coupling lies at the heart of a va¬ 
riety of novel electronic phen omena, including topologi¬ 
cal band insulators (TBI)PED an d Weyl semimetaliP. In 
a TBI the electronic band structure is topologically dis¬ 
tinct from that of a trivial insulator. The presence of zero 
energy edge modes, similar to the chiral boundary states 
of the quantum Hall effect, is a direct consequence of the 
topological order which defines a TBI. The predictiorP" 7 
and subsequent experimental verification 8 " 10 of novel 
topological phases in a variety of materials has spawned 
a vast new field of topological quantum matter. 

Recent developments in ultracold atomic gases al¬ 
low spin-orbit coupling to be tuned using Raman 
processes 11 13 l, expanding this frontier beyond electronic 
systems. Cold atom systems are already an excellent 
venue for the exploration of many-body physics 14 and the 
ability to simulate spin-orbit coupling promises to allow a 
detailed study of the interplay of spin-orbit couplin g an d 
strong interactions. Recent theoretical predictions^^ 
suggest that, on a lattice, the Bose-Hubbard model with 
spin-orbit coupling gives rise to a rich collection of effec¬ 
tive magnetic Hamiltonians in the strongly interacting 
limit that support a plethora of novel magnetic states, 
such as ferromagnets, antiferromagnets, spirals, and chi¬ 
ral textures. 

It is well known that systems with broken time-reversal 
symmetry exhibit off-diagonal Hall conduct ivitjf 17 !, and 
this transverse conductivity is nonzero in magnetic ma¬ 
terials even in the absence of an external magnetic field®. 
The study of the anomalous Hall effect has been ex¬ 
tended to a quantum geometrical viewpoint, with the 
Berry phase of the wavefunction taking a central role 
in this phenomena. The quantum anomalous Hall effect 
(QAH), a quantized version of the anomalous Hall effect, 
exhibits edge states which carry a quantized transverse 
conductivity, similar to the quantum Hall effect but with¬ 
out the requirement of an externally applied magnetic 
field. The origin of the QAH effect lies in an exchange 
interaction between a spinful itinerant particle and local¬ 
ized magnetic moments rather than from Landau levels 
due to the orbital effects of an external magnetic field. 
Although the QAH effect offers the promise of low dissi¬ 
pation transport without the need for an external mag¬ 


netic field, there are few experimentally accessible sys¬ 
tems which offer the combination of topologically non¬ 
trivial insulating behavior and magnetic structure. De¬ 
spite these challenges, the QAH effect has been observed 
in thin films of Cr-doped Bi 2 Te 3 , a magnetic topological 
insulatoPS. 

In this paper we propose a method for realizing the 
QAH effect in cold atom systems featuring itinerant par¬ 
ticles with spin-orbit coupling moving in a chiral spin tex¬ 
ture. We begin with a tight-binding model with Rashba 
spin-orbit coupling and a Hund’s coupling to a fixed mag¬ 
netic texture. After deriving an effective Hamiltonian 
for the system in the adiabatic approximation where the 
itinerant spin is always parallel to the local texture, we 
investigate some consequences of this effective Hamilto¬ 
nian for ferromagnetic, antiferromagnetic spiral and chi¬ 
ral textures. We show that chiral textures lead to a topo¬ 
logically nontrivial band structure for the itinerant parti¬ 
cles. By changing the strength of the spin-orbit coupling 
for the itinerant particle, we find the system undergoes 
a variety of topological phase transitions between vari¬ 
ous quantum anomalous Hall states. Finally, we propose 
that our results can be experimentally verified using Hall 
drift experiments and Berry phase interferometry. 

II. MODEL 
A. Hamiltonian 

We consider a single particle with two internal degrees 
of freedom, denoted f and J., on a two dimensional square 
lattice interacting with a localized spin texture through 
Hund’s coupling. We describe our system by the follow¬ 
ing Hamiltonian: 

-lAuii = t Si ’ Qa (Hi ) ab^'ib ■ 

{i,j) » 

(1) 

where d/t = (c^, ejjJ is a spinor of creation operators 
and the matrix = exp[iA • (r 4 — r,-)]. We see that 
on-diagonal elements of Rij describe spin-preserving hop¬ 
ping, while a non-Abelian gauge field A = (aa y , f3cr x , 0) 
leads to off-diagonal elements which describe spin flip 
hopping. We set /3 = —a in order to obtain the lattice 
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analog of Rashba spin-orbit coupling which we use for the 
remainder of this paper. We also take the lattice spac¬ 
ing to be unity such that for i,j which are neighbors, 
r, — Yj = 1 . The spin ai of the itinerant particle couples 
to the local spin texture denoted by with a coupling 
strength J > 0. 

We assume that the time scale of the itinerant parti¬ 
cle’s evolution is much shorter than that of the texture 
and so we consider the localized spins to be fixed in the 
spirit of the Born-Oppenheimer approximation. This can 
be accomplished in a cold atom system by creating the 
texture with a heavy species of strongly interacting par¬ 
ticles and allowing them to strongly couple with a lighter 
species of particles. The single particle approximation 
can be accomplished experimentally by taking the inter¬ 
actions between light particles to be exceedingly weak. 


B. Local Projection 



FIG. 1. A schematic illustration of the hopping processes 
which lead to the matrix elements t' and t". Due to the na¬ 
ture of the spin-orbit coupling, the hopping amplitudes now 
depend on the direction of the hopping as well as the orien¬ 
tation of the spin texture at each site. The faded red arrows 
show how a particle hops and then projects along the local 
texture and the dashed blue arrow represents the hopping. 
Particles hop without changing their spin and then project 
onto the local texture for t'. For t ", particles flip their spin 
from spin-orbit coupling while they hop to adjacent sites and 
then align with the local texture. 


In the limit of J/t —> oo, the coupling to the texture 
will split the spectrum into two bands and the low energy 
behavior of the system will be governed entirely by the 
lower band which is aligned locally with the spin texture. 
Because of this local alignment, we introduce operators 
representing the local orientation of the spin which 
creates a spin aligned with the local spin Sj at site i. 
In the usual manner, we can write these in terms of the 
creation and annhilation operators for spins aligned with 
the global ^-direction. We obtain the following unitary 
transformation 


ft 

•'it 


9cifij,(f>j) g s {9j , ( fij ) 

ft 





( 2 ) 


where g c (0j,(f>j) gives the component of the projection 
parallel to the global axis and g s (9j,(/)j ) gives the com¬ 
ponent of the projection antiparallel to the global axis. 
The orientation of the local spin at a site j is described 
by the polar and azimuthal angles (9j,(f>j) with respect 
to the global .c-axis. 

With Eq. ([2|, we can write the Hamiltonian in terms 
of the local spin creation and annihilation operators. For 
J/t —> oo, we can now project out the space where 
the spins are locally aligned with the texture by keep¬ 
ing only the terms in the Hamiltonian which are of the 
form 

The Hamiltonian now takes the form: 


//, 


proj 


e «,*+* cos(a) + t'i i+ s sin ( a ))fy i+ st ( 3 ) 


where S £ {±x, ±y} is a nearest neighbor of the site i. 
The hopping amplitudes f' . - and t”. - now encapsulate 
the information about the spin texture, the spin-orbit 
coupling and the direction S of the nearest neighbor. The 
explicit forms of the modified hopping matrix elements 


t' and t" are given in Appendix A. 

Embedded in this Hamiltonian is a real space vector 
potential which arises from the combination of the spin 
texture and Rashba spin-orbit coupling. 


C. Characterization of Spin Textures 

We fix a given spin texture a priori by choosing a set of 
spin configurations {0j, </>,}. In order to study the inter¬ 
play of spin-orbit coupling and coupling to a real-space 
spin texture, we investigate the effects of ferromagnetic, 
antiferromagnetic, spiral and chiral spin textures. Exam¬ 
ples of some of these textures are shown in Fig. Sard). 

These textures can be broadly characterized by their 
magnetic structure factor 

s(k) = ixy k - r *s j |, (4) 

3 

which is simply the Fourier transform of the local spin 
texture. In quantum materials it can be measured by 
neutron diffraction and in cold atom systems it can be 
accessed by Bragg spectroscopy®!. We show S( k) in Fig. 
[2|e-h) in the first Brillouin zone for a variety of textures 
on a 12 site by 12 site lattice. 

We also characterize a spin texture by its chirality. On 
a lattice, we define the local scalar chirality at a site i as 

Xi = (Si * (Sj-j-jz X S T Si * (Si_2i X S z— )) ’ (^) 

We will refer to \ = Xi as simply the chirality for 
a given texture. This quantity will obviously be zero for 
collinear or coplanar spin textures, such as ferromagnetic, 
antiferromagnetic and spiral textures. We call textures 
for which y / 0 as chiral textures. The chirality is a 
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topological invariant of the texture. 

Ferromagnet: We consider ferromagnetic (Fig. io 
textures oriented along the z axis. The magnetic struc¬ 
ture factor exhibits a peak at k = [0, 0] as shown in Fig. 

Hi. 

Antiferromagnet : In the antiferromagnetic phase (Fig. 
§>), spins point alternately along the ±z axis. The mag¬ 
netic structure factor has a maximum at k = [7T, 7r] as 
shown in Fig. (2f. 

Spiral: Spiral phases are composed of spins which cant 
with constant angle in one direction and which do not 
vary in an orthogonal direction. For simplicity, we con¬ 
sider spiral textures for which the spiral is oriented in 
either the x or y directions, rather than at an arbitrary 
angle. For a commensurate spiral oriented in the x ( y ) di¬ 
rection that winds with a period of L sites, the structure 
factor has a peak at [±27r/L, 0] ([0, ±27 t/L]) as shown 
in Fig. [2jx We call the angle ifj = 2n t/L between adja¬ 
cent spins the canting angle. We note that since spiral 
textures are coplanar, they have zero chirality. A spiral 
texture with a 4x1 unit cell is shown in Fig. [2^;. 

Meron: A non-coplanar texture in which the center 
spin in a unit cell is aligned with the z axis and spins at 
the edge of the unit cell take a polar angle of n/2 with 
respect to the z axis. Mapping the orientation of spins in 
a unit cell to their location on the unit sphere, one finds 
that the sphere is half covered. Merons have x = 1/2. A 
superlattice of 3x3 merons is shown in Fig. [2]l and its 
structure factor is shown in Fig. ©i- 

Skyrmion: Skyrmions are non-coplanar textures where 
the center spin in a unit cell is aligned with the z axis 
and spins at the edge of the unit cell are aligned with 
the — z axis. Mapping the orientation of spins in a unit 
cell to their location on the unit sphere, one finds the 
sphere to be fully covered. Skyrmions have X = 1 and 
they can be considered as a composite object composed 
of two merons. 


III. RESULTS 

We have exactly diagonalized the above single parti¬ 
cle Hamiltonian for several spin textures at a variety of 
spin-orbit couplings. After obtaining the single parti¬ 
cle spectrum, a wide variety of quantities can be inves¬ 
tigated. We derive effective Hamiltonians in the case of 
ferromagnetic, antiferromagnetic, spiral and meron tex¬ 
tures. We then characterize the topological nature of the 
bulk band structure by calculating the Berry curvature 
and the Chern number for each band. 


A. Projected Hamiltonian in Simple Textures 

For a purely ferromagnetic texture, t ". is zero for all 
i and t' - is unity for all i. It is clear that by projecting 
into the locally aligned basis, systems with ferromagnetic 


textures are similar to a free particle for closed boundary 
conditions or the familiar particle in a box for those with 
open boundary conditions. Spin-orbit coupling leads to 
a modulated effective hopping strength 

A FM (k) = —fcos(a)[cos(fc x ) ±cos(fc y )] (6) 
shown in the left panel of Fig. [3] 

Due to the opposite orientation of neighboring spins in 
the antiferromagnetic texture, it is clear that t[ • is identi¬ 
cally zero for all values of spin-orbit coupling a. Without 
spin-orbit coupling, this system is trivial, since the local 
projection destroys hopping amplitudes between oppo¬ 
sitely aligned spins. The spin-flip terms in the Hamilto¬ 
nian in Eq. ([I]) dominate the contribution to the hopping 
elements in projected hopping elements in Eq. © and 
lead to an effective Hamiltonian for the antiferromagnetic 
texture 

Haf = -2t sin(a) $]^AF(k)$ kt . (7) 

k 

Here J^b) where the /^’s are defined as 

the Fourier transform of the operators defined in Eq. © 
and the A or B refers to the two sublattices of the anti¬ 
ferromagnetic texture. The kernel %AF(k) in Eq. (|7|) is 
given by 

%AF(k) = sin(A; x )T x ± sm(k y )T y (8) 

where r are the pseudo spin Pauli matrices associated 
with the A and B sublattices. We find the energy disper¬ 
sion shown in the right panel of Fig. [3] to be 

E^ f ( k) = ±2i sin(a)^/sin 2 (fc x ) +sin 2 (k y ). (9) 

This spectrum is identical to that of the Hamiltonian for 
a free particle with Rashba spin-orbit coupling and no 
spin-preserving hopping. Here we find Dirac cones at the 
time-reversal invariant momenta. We note that since we 
do not obtain a gapped band structure, we cannot define 
a topological invariant for the antiferromagnetic case. 

For a spiral texture, we pick the direction of the spiral 
canting to be in the x-direction without loss of general¬ 
ity. Motivated by recent theoretical predictionf®, we 
investigate specifically the case of L = 4 spiral. For 
this texture, the band structure does not have a sim¬ 
ple closed analytic form. Notably, both t' and t" are 
non-zero, leading to a nontrivial dependence of the dis¬ 
persion on a. In Fig. [4j we show the band structures 
for a = 0, 7t/4 and n/2 across the reduced Brillouin zone 
{—7T < k y < 7T, —7t/4 < k x < 7t/4}. We draw particular 
attention to the flat bands at zero energy for a = 7r/4 
where commensuration of spin-orbit coupling and the an¬ 
gle of the spiral lead to destructive interference which 
eliminates hopping along the direction of the spiral. 
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FIG. 2. (a-d) The real space spin configurations for several spin textures: (a) a ferromagnetic texture, (b) an antiferromagnetic 

texture, (c) 4x1 spirals on a 12x12 lattice, (d) 3x3 merons on a 12x12 lattice, (f-h) The spin structure factors for several textures 
on a 12 by 12 site lattice: (e) a ferromagnetic texture, (f) an antiferromagnetic texture, (g) spirals with a 4x1 unit cell in the 
y direction, (h) sixteen 3x3 merons 



FIG. 3. (Left) The band structure for a particle in a ferromagnetic texture is identical to that of a free particle with a 
cos(a) modulated bandwidth. (Right) The bands for an antiferromagnetic texture are identical to that of the Rashba portion 
of the Hamiltonian in Eq. |l]) with spin-preserving hopping with a sin(a). There are Dirac cones at the time-reversal invariant 
momenta for the antiferromagnetic bands. 


B. Chern Numbers and Berry Curvature 


where the Berry connection A„ (k) is defined in terms of 
the wavefunction of the n th band |n(k)) aiP^ 


The band structure of a periodic system determined 
by a Bloch Hamiltonian iT(k) can be characterized by 
the Chern numbeJ^tlSl. The Chern number for the n th 
Bloch band is defined as 

c n = / rf2k ’ Vk x A «( k ) 


A n (k) = (n(k)| *V k |ra(k)). (11) 

The Hall conductivity is given by 

e 2 x - 

a xy = / , Cw; (12) 


( 10 ) 


n 
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FIG. 4. The band structure in the reduced Brillouin zone for a 4x1 spiral texture with a = 0 (left), a = 7t/4 (middle) and 
a = 7t/2 (right). We note the presence of zero energy bands for a = tt/4. 



where the summation index n runs over the filled bands. 

The presence of a spin texture in the Hamiltonian in 
Eq. 0 explicitly breaks time-reversal symmetry in gen¬ 
eral and so the Chern number c n of the n th band of H pm j 
in Eq. © will not be zero in general. A square spin 
texture with sides of length L in a system with periodic 
boundary conditions is described by a Bloch Hamiltonian 
which can be represented as an L 2 dimensional matrix in 
momentum space. Due to the numerical nature of our ex¬ 
act diagonalization, we can diagonalize our system only 
on a discrete mesh of points ki within the Brillouin zone. 

Following Fukui, Hatsugai and SuzukP^, we introduce 
a U(l) link variable 

U n ,ii{ki) = (n(ki)\n(ki + p,))/\{n(ki)\n(ki + /t))| (13) 

and a lattice field F analogous to the continuum Berry 
curvature Vk x A n (k): 


F n (ki) = In 


/ Un^x^kl^Un^y^ki T x) \ 

\Un,x{ki + y)U n , y (ki) J 


(14) 


where we choose F(ki) to lie within the principal branch 
of the logarithm. In the logarithm in F(ki), the U(l) 
gauge field is summed around a single plaquette starting 
at the momentum k/. Lastly, we take the lattice Chern 
number to be 


= ( 15 ) 

i 


where l — 1 , ...,L 2 . In addition to providing a robust 
method of calculating the Chern number for a modestly 
sized mesh of points in momentum space, the lattice field 
F n is manifestly gauge invariant. 

In Fig. 0, we show the sum of the Chern numbers for 
the lowest two bands 


ri 2 = ci + c 2 (16) 

as a function of spin-orbit coupling calculated for a parti¬ 
cle in a 3x3 meron texture. In Fig. 0b,c,d), we show the 
energy dispersion E( k) along the principal axes of sym¬ 
metry in the Brillouin zone for various values of spin- 
orbit coupling. We see that tuning the spin-orbit cou¬ 


pling causes the Chern number n 2 to change, signaling 
a topological phase transition. Each time n 2 changes 
the gap between the second band and the third band 
closes. For a range of a there exists a global gap be¬ 
tween the second and third bands and in such cases the 
chemical potential can be placed within the global band 
gap giving n 2 that is proportional to the Hall conductiv¬ 
ity cr X y Recent experiments have shown that the Chern 
number can be measured in cold atomic gases 26 by using 
the transverse deflection of an atomic cloud in reponse 
to an optical dipole force. 

Rather than considering only the Chern number, which 
gives us information about the presence of a band cross¬ 
ing, we can consider the monrentunr-resolved quantity 
F n across the Brillouin zone. Peaks in this quan¬ 
tity reveal the locations of crossings of bands in mo¬ 
mentum space. These can be experimentally detected 
by Aharonov-Bohm interferometry in momentum space. 
Recent experiments 2 ' have shown that the Berry phase of 
two Bose Einstein condensates traversing separate paths 
in the Brillouin zone can be measured with high momen¬ 
tum resolution. The topological phase transitions that 
we predict can be measured by performing momentum- 
space Aharonov-Bohm interferometry of two BECs sent 
along paths in the Brillouin zone which enclose the asso¬ 
ciated Berry curvature. 

In Fig. 0e,f,g), we show the lattice Berry curvature F 2 
of the second band for several values of spin-orbit cou¬ 
pling. For a = 0 (Fig0), we see one a very large negative 
contribution to the Berry curvature at k = (7 r, 7r), arising 
from an isolated band crossing that acts like a monopole 
of Berry curvature. For a = 0.167T (Fig0), rather than 
a monopole, we see a ring of Berry curvature through 
a portion of the Brillouin zone, indicating a band cross¬ 
ing that occurs for a locus of momenta points. When 
a = 7t/2 (Fig [5^) we see a positive peak in the Berry 
curvature for k = ( 0 , 0 ) and four other near the corners 
of the Brillouin zone. 


IV. CONCLUSIONS AND DISCUSSION 

We have shown that the properties of particles moving 
in spin textures can be further modified by spin-orbit cou- 
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FIG. 5. (a) Sum of Chern numbers n 2 as a function of a for the lowest two bands of a system with unit cells of 3 by 3 merons. 

These lowest two bands are well separated from the others except for 0.5 < a < 0.67 and 1.2 < a < 1.4 and so, excluding 
these regions, the Chern number can be interpreted as the Hall conductivity. Spin-orbit coupling has a non-trivial effect on 
the topological nature of the bands for a fixed texture. (b,d,c) Band structure for 3 by 3 merons, along the principal axes of 
symmetry T —> X —> M —> I\ (b) a = 0, (c) a = 0.167T, (d) a = 0.57T. (e,f,g) Berry Curvature F 2 for the second band of a 
3 by 3 meron system for (e)a = 0, (f) a = 0.16-7T and (g) a = 7t/ 2. Peaks in the Berry curvature identify values of k which 
contribute to the Chern number. 


pling in several important ways. In simple textures, such 
as ferromagnets and antiferromagnets, the spin-orbit cou¬ 
pling can be used to control the bandwidth of the particle 
while keeping the other properties of the band structure 
invariant. We have shown that for a 3x3 meron tex¬ 
ture the Chern number is nontrivial for the lowest two 
bands. Above these two bands either a local gap or a 
global gap exists for nearly all values of a signaling a 
wide range of Chern metal and Chern insulator states. 
By simply changing the value of spin-orbit coupling of 
the itinerant particle, the Chern number gets modified, 
signaling transitions between topologically distinct quan¬ 
tum anomalous Hall states. 

Recent work in quantum materials has focused on the 
current driven motion of skyrmion^D and on the emer¬ 
gence of a finite anomalous Hall conductivity of s- elec¬ 
trons coupled to chiral magnetic textures of localized 
d- electrons 29 . We emphasize that our model predicts 
a quantum anomalous Hall effect whose Chern number 
is tunable by simply changing the strength of Rashba 
spin-orbit coupling of the itinerant particle, a qu antity 
which is very accessible in cold atom experiments 13 30 . 


We have shown that multiple signatures of these topo¬ 
logical phase transitions are experimentally accessible 
via Bohm-Aharonov interferometry and Hall drift exper¬ 
iments. 

Rather than fixing a texture a priori, in future calcu¬ 
lations we will use Monte Carlo simulations to generate 
thermally fluctuating magnetic textures to explore the 
effect of temperature on the Berry curvature and the 
robustness of the Chern number. Another aspect that 
could be important is the influence of the kinetic motion 
of the itinerant particle on the magnetic texture itself. 
Cold atoms also provide an exceptional setting for the 
study of strongly correlated particles and the associated 
fractional quantum Hall physics on a lattice without the 
inclusion of an external magnetic field. We have demon¬ 
strated that tuning between various quantum anomalous 
Hall states is possible by varying the spin-orbit coupling 
of the itinerant particle and it is possible that this tun- 
ability will allow for a controlled way of exploring the 
transitions between various fractional quantum anoma¬ 
lous Hall states. 

We thank W. S. Cole for useful discussions. T.M.M. 
and N.T. acknowledge funding from nsf-dmrl309461 
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V. APPENDIX A: DETAILS OF THE 
PROJECTED HOPPING ELEMENTS 


J/t —> oo, the coupling to the texture will split the spec¬ 
trum into two bands and the low energy behavior of the 
system will be governed entirely by the lower band which 
is aligned locally with the spin texture. Because of this lo¬ 
cal alignment, we introduce operators representing the lo¬ 
cal orientation of the spin which creates a spin aligned 
with the local spin at site i. In the usual manner, we 
can write these in terms of the creation and annhilation 
operators for spins aligned with the global ^-direction. 
We obtain the following unitary transformation 


ft 

Jjj 


9c{Qj i ) 9s(@ji <t>j ) 

ft 

. Mi . 


gliOj.Oj) -g* c {9j,(f)j) _ 


Mi 


(17) 


where g c (9j,cj)j) gives the component of the projection 
parallel to the global axis and g s (9j,(j)j) gives the com¬ 
ponent of the projection antiparallel to the global axis. 
The orientation of the local spin at a site j is described 
by the polar angle 9j and azimuthal angle <f>j with re¬ 
spect to the global z-axis. In Eq. the elements of 

the rotation matrix are given by 


gc{0j,<t>j) = cos(6t,/2)e 

(18) 

9s{0jAj) = sin(0j/2)e*^ /2 . 

(19) 


In order to derive the projected hopping Hamiltonian in 
Eq. @, we rewrite all operators in Eq. 0 in terms 
of the local operators given in Eq. |L7| ) and keep only 

terms of the form fj fj. We find that the term t' ? 

i,i+S 

originates from the hopping without spin-orbit coupling 
in the global frame and is found to be 


t'i,j = 9c(0i, (t>i)g*(9j, (f>j) + g s {9 i ,(l)i)g*(9j,(j) j ) 

= cos( 0 4 /2) cos(^/2)e- iWi -^ )/2 

+ sin(0,;/2) sin(0 i /2)e iWi -^ )/2 (20) 

for any neighbor j. The spin-orbit coupling in the global 
frame results in the other hopping amplitude t" . - which 

is found to be 


ti,j=i+y = - i (9s(S l ,c/) i )g*{9 j ,(l)j)+g c (9i,(l}i)g*{9 j ,(l) j )) 

= — *(sin(#i/2) cos(0,-/2)e*^ i+ ^.i)/ 2 

+ cos(6t/2) sin(6» J /2)e- i( ^+^ )/2 ) (21) 


for hops up, 

= i{9s{0i, ^g^Oj,^) + gdOiA^gt^j^j)) 

= z( sin( 0 i/ 2 ) cos( 9 j 

+ cos(0j/2) sin(%/2)e- i W <+ ^ )/2 ) (22) 


In this appendix, we show the explicit form of the ef¬ 
fective hopping amplitudes t! and t". In the limit of 








for hops down, 


for hops left. 


t",j=i+x = OsiOiA^gc^j^j) - 9c(Gi,<Pi)g* s (0 j ,(f> j ) 

= sin(0.j/2) cos(0j/2)e i ^ i+ ^^ 2 

- cos(6» i /2) sin(0 i /2)e- i( ^ + ^ )/2 (23) 


for hops right, and 


t",j=i-x = -(gsiOi^iWciOjAj) ~ 9c{Qu<!>i)9* a {6j,<t>jj) 
= cos(0,;/2) sin(0j/2)e _ *^ i+l ^/ 2 

- sin(0,/2) cos(0 i /2)e iW * + ^ )/2 (24) 



